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Abstract 

We derive the kinematical constraints which characterize the decay of any massless 
particle in flat spacetime. We show that in perturbation theory the decay probabilities of 
photons and Yang-Mills bosons vanish to all orders; the decay probability of the graviton 
vanishes to one- loop order for graviton loops and to all orders for matter loops. A general 
power counting argument indicates in which conditions a decay of a massless particle 
could be possible: the lagrangian should contain a self-coupling without derivatives and 
with a coupling constant of positive mass dimension. 
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The massless particle which we best know, the photon, is certainly stable for very long 
periods. The experimental evidence concerning the properties of the neutrino (admitted it is 
really massless) is less strong, but it is generally regarded as stable too. 

Nevertheless, kinematics allows in principle the decay of a massless particle, provided the 
products are massless and their momenta have the same direction and versus of the initial 
momentum (compare Section 1). This means that the Mandelstam variables of the process 
vanish, so that its amplitude, regarded as a function of Mandelstam variables, must be computed 
in this particular limit 0. Moreover, even if the limit of the amplitude is not zero, the phase 
space for the products reduces to a line in momentum space and therefore its volume tends to 
vanish. 

An almost exact "collinearity" of the products is usually observed in the decay of any 
particle for which E ^ m. Four- momentum conservation implies (see for instance fl]]) that 
the mass of a particle produced in an annihilation process is proportional to the sine of the 
angle between the momenta of the colliding particles; conversely, in the limit m — the decay 
produces coUinear particles. 

Limits of this kind (m ^ 0) are common in the treatment of infrared singularities in quan- 
tum field theory (f^; see also Sect. 5). In our case, however, the assignment of an infinitesimal 
mass m to the particles involved in the decay is unsuitable as a regularization technique. In 
fact, let us consider the decay of a massless particle of energy E into n coUinear particles of the 
same kind, with energies Ei such that J2i=i...n Ei = E. This process is kinematically allowed 
(compare Section 1; n must be odd if the initial particle has nonzero helicity), but if we give 
the particles an infinitesimal mass, it becomes obviously impossible (suppose to observe it in 
the rest system of the initial particle). 

We shall then work from the beginning with massless particles and introduce a different 
regularization, involving a weak external source J which gives the initial particle an infinitesimal 
additional energy (and/or momentum) uj, so that its 4-momentum is put slightly off-shell. This 
regularization technique proves to be quite effective, as it also allows an estimate of the decay 
probability by power counting. 

Let us now come to the specific cases we treated. In QED it is possible to show in a general 
way by means of the Ward identities that the decay amplitude for 7 — 71 + ... + 7„ {n odd) is 

^In the four-particle amplitude we mean by Mandelstam variables the usual ones, s, t, u; for amplitudes 
with more external massless particles, they are taken to be all the possible scalar products between the external 
four- moment a . 
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a symmetrized sum of terms which can be factorized into a finite scalar part and a tensor part 
that vanishes when all the external momenta are aligned. An analogous reasoning holds for the 
neutrino. In both cases, it is crucial that the loop amplitudes contain in the denominator the 
masses of the fermions or of the vector bosons, respectively. 

Another example of massless particle is the graviton. Here we do not have any experimental 
evidence yet. It has been suggested [Q that the non-linearity of Einstein equations could lead 
to a "frequency degeneration" in gravitational waves, a phenomenon which from the quantum 
point of view would correspond to a decay of the graviton into more gravitons of smaller energy. 
We were able however to prove through a generalization of the procedure applied to QED that 
the amplitude of this process vanishes in perturbation theory around the fiat background. In 
this case the negative mass dimensionality of the Newton constant plays a role analogous to the 
fermion masses in QED. At the non perturbative level, the hypothesized existence of a small 
scale cosmological constant could change the situation (see below). 

The case of the gluon, although physically quite academic due to the confinement, is par- 
ticularly interesting because the amplitude of the decay g gi + ... + gn {n odd) is finite for 
n = 3 and divergent for n > 5. (The Ward identities still allow a factorization of this amplitude 
like in QED, but the scalar parts now contain poles.) Nevertheless, the total decay probability 
is zero because the phase space for the products is suppressed strongly enough to compensate 
for the divergence in the amplitude. We thus have here a typical example of cancellation of 
infrared divergences in the computation of a physical quantity. 

A general power counting argument indicates in which conditions a real decay of a massless 
particle could be possible: the lagrangian should contain a self-coupling without derivatives 
and with a coupling constant of positive mass dimension. This is precisely what happens in 
quantum gravity in the presence of a cosmological constant, and in fact it has been suggested 
that in this theory strong infrared effects could become relevant But one must remind that 
in the lagrangian the cosmological constant also multiplies a term which is quadratic in the 
field and thus generates an effective mass for the graviton (if A < 0) or an unstable theory (if 
A > 0) p. A possible way to elude the problem is to admit, like in lattice theory, that the 
effective cosmological constant vanishes on large scales but not on small scales and is negative 
in sign (compare Section 5). This latter approach is however out of the scope of our paper. 

The structure of the article is the following. Section 1 is concerned with kinematics. In 
Section 1.1 we give a list of simple kinematical properties which characterize the decay of any 
massless particle. These properties are only due to Lorentz invariance and to the conservation 
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of the total four-momentum and angular momentum. In Section 1.2 we reexpress in a more 
manageable form the Lorentz-invariant decay measure defined on the phase space of n massless 
product particles, under the condition that also the initial particle is massless; specializing to 
the case n = 2 we compute explicitly the lowest order decay probability in the toy-model scalar 
theory. In Section 1.3 we introduce an infrared regularization which allows the computation 
of the decay amplitudes in the limit of vanishing Mandelstam variables. In Section 2 we give 
a dimensional estimate of the decay probability of the photon, the neutrino, the gluon and 
the graviton. After recalling in Section 3 how the exact proper vertices are connected to 
the complete perturbative expression for the decay amplitude, in Section 4 we use the Ward 
identities for QED, Yang-Mills theory (YM) and Einstein quantum gravity (QG) to give an 
estimate of the regularized amplitudes. In Section 5 we comment on the relation between the 
infrared singularities which occur in our computations and the usual infrared singularities of 
quantum field theory. Finally we present a few brief speculations about the possible role of a 
non-vanishing cosmological constant in the decay of the graviton. 

1 General kinematic properties. 
1.1 Consequences of Lorentz invariance. 

We list here the most general properties of the decay of a massless particle. They are due only 
to the Lorentz invariance of the process and to the conservations of the total four-momentum 
and angular momentum. As we mentioned in the introduction, some of them can be proven 
taking the limit m ^ in the corresponding formulas for massive particles |jl|. Properties 1, 
2, 3, 6 can also be found in ref. 0. 

Property 1. - A massless particle can only decay into massless particles. - In fact, through 
a suitable Lorentz boost we can make the energy of the initial state arbitrarily small. If, per 
absurdum, in the final state massive particles were present, the energy of this state would be in 
any reference frame equal or bigger than the sum of the masses. 

Property 2. - Let us suppose that the impulse p ° of the initial particle is oriented in a certain 
direction and versus, for instance let its four-momentum have the form 

/ = E', 0, 0) (1) 

Then also the impulses p^...p" of the n product particles are oriented in the same direction 
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and versus; in our example we shall have (Fig. 1) 



n 



= {E\ E\ 0, 0); 



i = 1, ...,n; 



J2e' = e'>. 



(2) 



In an arbitrary Lorentz frame this can be rewritten as 



(3) 



where 1 > Aj > 0, i = 1, 2, n, J2i K = 1- 



E^ 





■> 



E 



Figure 1: CoUinearity property (Property 2). 



- Also this property depends on the fact that through a suitable Lorentz boost along z we 
can make the energy of the initial state arbitrarily small; while if per ahsurdum in the final 
state some transversal momenta were present, their contribution to the energy would not be 
affected by the boost. 

Property 3. - If the initial particle has helicity h and decays into n particles of the same 
helicity, n must be odd. - The proof follows directly from Property 2 and from the conservation 
of the angular momentum. 

Property 4. - In the decay of a massless particle, all the scalar products {p'^-p'), i,j = 0, l,...,n 
vanish. This means that the Mandelstam variables vanish. - The proof follows directly from 
Property 2. 

Property 5. - If represents the polarization vector of the i-th particle involved in the decay, 
in a gauge such that (p* ■ e*) = 0, then we have also (p* ■ e^) = for i,j = 0, 1, ...,n. - Once 
more, the proof follows directly from Property 2. 

Property 6. - If a massless particle decays, its lifetime r in a reference frame where its energy 
is E^ has the form 



where ^ is a constant which depends on the dynamics of the process and has dimension [mass]^'^. 
- This property holds also for massive particles, for which the constant takes the form ^ = 
Trest/^- The proof is elementary (see for instance 0). 



lO 



(4) 
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1.2 The decay phase space measure d^n' 

We recall that according to quantum field theory the decay probability (per unit time) should 
be computed by the general formula 



T 



r 



2^0 



n>2 



En 



{2'Kf2E 



(5) 



where T„ is the quantum amplitude for the decay process into n product particles of momenta 
{p*}. If the final particles have helicity or internal quantum symmetry numbers T„ includes the 
sum over these degrees of freedom. 

Actually, both eq.s and give physically realistic predictions as far as: 

(1) the energy uncertainty /\E of the first particle fulfils the condition /\E ^ E] 

(2) [] the finite energy resolution e of the decay detector can be neglected. In general the 
detector will be unable to recognize a decay process in which one of the outcoming particles 
has energy E' such that E — E' <^ e. In order to compute the correct detection probability 

one should in principle subtract from formula the total probability of all events of this 
kind. Nevertheless, for the theories considered in this paper one finds that this effect is indeed 
negligible (in perturbative QED, YM and QG we will find F = 0, whence it follows = 0, 
since F^ < F; in the toy-model at order O(A^) considered below one finds that F — F^ ~ e). 
We plan to devote more attention to the general issue elsewhere, by considering examples of 
theories for which condition (2) is not fulfilled. This requires an approach to IR divergences as 
in the Kinoshita-Lee-Nauenberg theorem 0. 

A closer look at the measure appearing in the integrals on the RHS of formula (^ is now 
very useful. When all particles are massless, it is possible to express the Lorentz-invariant decay 
measure 



dij'n = n 



in the following form: 



2a 



2n-2 



E^ 



\ {2'Kf2E'' 



X{d'p'5\p^)e{pi] 

.i=l 



(6) 



1=1 



5 [p' -T.pl 



2 



i=l 



1 



Ept 



\i=l 



i=l 
n-2 



(7) 



^We thank M. Abud for drawing our attention to this point and to several further subtleties which are 
required by a correct physical interpretation of infrared divergences related to massless particles 
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where pi and denote the longitudinal and trasversal part of the impulse p * with respect 
to the direction and versus identified by p °, 9 is the step function, and the adimensional 
coefficients a are those which appear in the expression of 5^{x) in polar coordinates: 



=a-^5(|f|')|f| 



2-m 



The 5-functions occurring in formula (J^) show that the support of dn^ is concentrated around 
(the infinitesimal neighbourhood of) the collinearity region, which is characterized by all sets 
{p*} satisfying relation (H). 

The collinearity property (|]) follows from the sole condition p° — Yh=i = [imposed by 
the (5-function contained in formula (^] if all are null vectors. In fact, we observe that 



\i=l 



i=l 



= [T./) =(p°)' = o- (9) 

The 3- vector Z^ILiP * has length i < J2i=i |p and the equality holds only if p * = Xip, for 
some p and some {Aj} all of the same sign; inserting this into the relation p° — J27=iP^ = we 
find eq. {§. 

The squared amplitude |T„p depends only on the Lorentz invariants (p* ■ p^)- But in the 
collinearity region p* ■p-' = 0. Thus in this region a finite |T„p may only be a function (possibly 
trivial) of the invariants Aj defined in formula (|^); in this case the corresponding integral can 
be easily performed and gives a finite (possibly vanishing) result. In particular, when n > 3 the 
factor in the square bracket of formula (|^) is set equals to zero by the 5-functions, and thus if 
the amplitude of the decay is finite, the corresponding total probability is zero. If |T^| diverges, 
we may introduce a suitable regularization in order to make the integration easier (see Section 
1.3). 

For a massless scalar field theory with self-coupling of the form A0^ the phase space integral 
(^ with n = 2 coincides, up to a factor A^, with the probabihty of the decay of a particle into two 
particles, computed perturbatively to lowest order. This is a concrete example of computation 
of a finite decay probability, although with the known limitations of the X(f)^ theory fl 

Setting = 2 in d^) and performing the integral (the square amplitude does not depend on 
p^, p^ and is equal to A^) we obtain 

dl^2 = ^J dp'dp' 5[E' - (p'Y - ip'Y]9[ip'Y]9[ip'Y] = 2a,. (10) 



"it is known that the action is not hmited from below and that the radiative corrections do not preserve 
TO = 0. 
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The present conclusion that F is finite to order coincides with that of the dimensional 
analysis (considered in Section ||) applied to this case, in which the coupling constant has 
positive mass dimension. 



1.3 Regular izat ion through an external source. 

We would like now to introduce an infrared regularization in order to allow a quick estimate 



coUinearity region). 

Obtaining such a regularization is not trivial. The most common infrared regularization 
technique, which consists in giving the soft particles a small mass ^ which eventually goes to 
zero, does not work in the present case, because the (regularized) process in which one particle 
of mass fi decays into more particles of the same mass has obviously zero probability. 



Instead, a better approach is to put external momenta slightly off-shell in a way controlled 
by an infinitesimal parameter u. 

Let us suppose (Fig. 2) that a very weak external source J gives the decaying particle (state 
I) an infinitesimal additional energy u. The exact nature of the source and of the particle which 
carries the energy u are not essential. For instance, if J represents a classical field, the energy 
can be carried by an on-shell boson with four-momentum [uj, 0, 0, uj); by absorbing the boson, 
the initial particle gains a small transversal impulse (state II) too. Alternatively, the energy 
UJ could be carried by an off-shell boson produced in J through an annihilation process, with 
four-momentum {u, 0, 0, 0); by absorbing the boson, the initial particle gets off shell too. More 
generally, we will assume that after the interaction the four-momentum of the initial particle 
will have the form 



of the integrals (^ in all cases (including the case in which the amplitude T„ diverges on the 





Figure 2: Factorization of the decay amplitude. 



p f = p' 




(11) 
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where W is adimensional and (p")^ = 0. It is not necessary to make any special assumption on 
the four- vector 6° at this stage (in Sections 3, 4 we will prefer to specialize the discussion by 
assuming that Ip ■ = 0, and other similar conditions for the p*'s). 

At this point the decay takes place; the products (state III) have now a small tranversal 
impulse of order uj and the Mandelstam variables (p* -jP) are of order uj (at least). The partial 
decay probability into n product particles is written as a sum over intermediate states (compare 
(eq. I) 

where T is the appropriate evolution operator. When u; ^ 0, the factor tends to 

which is the dependence that we expect on the basis of Lorentz invariance (compare Property 
6). Thus in this limit the integral In appearing in the preceding formula does not depend on 
E. Summing up we obtain 

the only the only massive parameters on which /„(co') depends are uj and the massive parameters 
possibly present in the theory that we are considering. This allows in most cases to estimate 
dimensionally whether r„ is finite, vanishes or diverges in the limit uo Q [note that the mass 
dimensions of /„, \{II^^^\T\III)\^ and of djjin are respectively equal to 2, 2(3 — n), 2n — 4]. We 
shall give some examples of this in the next Section. 



2 Power counting. 

In several cases the integral /„ can be estimated by simple arguments (often dimensional con- 
siderations alone are enough). 

For instance, in QED the four-photons amplitude is given to lowest order by the four 
fermions loop (fig. 3a). It is easy to realize that the loop integral gives a 4-th degree homoge- 
neous polynomial in the dimensionless variables ^ 0, where m/ is the mass of the fermion. 
The integral Is will therefore be proportional to 

where a is the fine structure constant. All behaves as though 

T^^a'(—\ . (15) 
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To be precise, the behaviour T3 ~ u'^ holds only for some specific choices of the "slightly off- 
shell" external momenta p*, whereas in any case T3 = 0(ct;^) at least; the integration transforms 
the remaining dependence of T3 on p*, if any, into an additional uj"^ factor. 

The above result can be generalized to the n-fermions loop: the key point is that the 
fermionic propagators of the loop produce masses in the denominator. The case of the neutrino 
is analogous: the masses of or appear at the denominator in the amplitude. In both 
cases, since the amplitude is proportional to a positive power of the regularizator u, it vanishes 
in the infrared limit due to (|12|). 




Figure 3: (a) Fermions square loop, (b) Gravitons or gluons loop. 

In the case of pure quantum gravity we have tree and one-loop graviton diagrams with k 
external legs (fig. 3b). Explicit expressions for the k = 4 amplitudes have been given by H. 
In any case, these amplitudes contain positive powers of the constant k = ^/l6^^G and then, 
like in QED, they behave always like a positive power of u and cause the decay probability to 
vanish. 

In the case of QCD the amplitudes do not contain dimensional constants. We expect that 
the decay amplitude of the gluon into three gluons, being adimensional, tends to a constant 
when — > 0, and this is in fact what happens 0. The decay amplitudes of a gluon into 5, 7 
... gluons have mass dimensions -2, -4 ... respectively, so they diverge when — 0; but this 
divergence is compensated in the phase space integral by a bigger positive power of u in such 
a way that the probability behaves like cu'^/E^^^ and thus vanishes in the limit. 

We are not going to apply this power counting argument to all possible theories and cou- 
plings, since it is in each case quite immediate. As a last example, we may wonder whether 
a photon can in principle decay due to the gravitational interaction, through diagrams with 
external photons and one loop of gravitons. Since the coupling constant k has mass dimension 
-1, while the fine structure constant a is adimensional and there are no masses involved, we 
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conclude once more that the amphtude of the process vanishes in the infrared hmit. 

It is clear from the discussion above that a r„ different from zero can be only obtained 
when the square amplitude is proportional to a sufficiently high negative power of to. Since in 
perturbation theory the coupling constants always appear in the numerator, this means that 
the amplitude must contain a coupling constant with positive mass dimension. We shall return 
on this point in the conclusions. 

3 Diagrammatics: cc;-dependence of the decay amplitudes. 

The dimensional arguments of the previous section determine the a;-dependence of the decay 
probability only for the pure gauge theories (YM, QG), where the only parameter in the action is 
the coupling constant. If additional dimensionful parameters appear in the action (as it happens 
for instance when the gauge field is coupled to some massive field) the previous arguments, as 
we have seen in the QED example, must be completed by some additional information. In 
general, a more explicit analysis of the perturbative expansion and use of Feynman diagrams is 
therefore needed in order to estimate the total decay probability. In this and in the following 
section we carry it out in such a way to determine not only the a;-dependence of the total 
decay probabihty, but also of the decay amplitudes (i.e., of the probabihties of the single decay 
channels) . The general results for the former will be essentially the same as those found by the 
dimensional arguments in section 2. Thus, we conclude that the decay probability of the gauge 
bosons of QED, YM, QG vanish. 

Before starting, let us define a "decay configuration" as follows: it is a pair of {n + 1) four- 
momenta and (n + 1) polarization vectors (p*, £*)i=o,i,...,n satisfying the properties {piY — 0, 
Yh=qP^ = 0, (e* ■ p") = 0, Pq > 0, Po < for / = 1, ...,n. We thus agree that the signs of 
the four-momenta of the outgoing particles are reversed. As we have seen, for particles with 
non-zero helicity n must be odd. 

We will start the analysis of the perturbative expansion from the tree level: a sum of 
truncated connected tree-diagrams with {n-\-l) external lines will give the lowest order (in 
h) contribution to the decay amplitude of 1 gauge boson in n gauge bosons. Higher order 
corrections will involve truncated connected diagrams with one or more loops. To formally 
compute the "exact" decay amplitude one has to replace in each tree diagram every boson 
propagator with the corresponding exact boson propagator, and each m-boson vertex with the 
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corresponding m-boson proper vertex (i.e. one-particle-irreducible Green function)]^. In order 
to get the /i^-order approximation of the decay amplitude, one simply has to retain the terms 
of order < r in this formal " exact" expression. As we will see, the Ward identities imply 
that when approaching a decay configuration: (1) in QED the decay amplitude of a process 
with m external photons vanishes; (2) in QG the decay amplitude of a process with m external 
gravitons or photons vanishes; (3) the decay amplitudes of processes with external Y.M. bosons 
may be finite or diverge, but in such a way that the corresponding decay probabilities vanish. 

3.1 Tree level 

Let us start from the Feynman vertices with m gauge massless bosons (m > 3) [see the actions 
(^51)]: we draw them in fig. @. The diagrams are to be understood as truncated in the external 
lines. In QED there is no m-photon vertex. In YM there are only two m-gluon vertices (for 
m = 3,4). In pure QG there is one m-graviton vertex for every m > 3; if coupling of gravity 
with the electromagnetic or the Yang- Mills fields is considered, then there are also vertices with 
k spin-1 bosons (photons or gluons) and r gravitons, for A; = 2,3,4 and r > 1. In the figures, 
a wavy line in the QG case will denote either a graviton or another gauge boson (a photon or 
a gluon). 




Figure 4: Feynman vertices 

At the tree level, the decay amplitudes '^^ T^^^ Tf ... of YM, QG are the sum of 
the diagrams in fig. (^. 

Tree diagrams involving ghost lines do not contribute to T^^*^*^. In fact, even though ghosts 
are massless, diagrams with external ghosts are zero when multiplied by physical polarization 
**In principle, propagators and proper vertices could be computed even in two different gauges, in order to 
simplify calculations, see Ref. jl^ 
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Figure 5: Tree level amplitudes: (QG) means that the diagram in Tg'''^^ is present only in QG. 

vectors, and diagrams with internal ghost lines (propagators) have necessarily also external 
ghost lines, by ghost number conservation. One can easily verify that in QG the decay ampli- 
tudes with only m external gravitons or photons vanish {T^^^ = 0) in any decay configuration, 
because each vertex is quadratic in the momenta fc*, implying an overall (fc)^ dependence of 
each separate diagram in fig. @; when contracted with the external polarization vectors, this 
will give zero, since in the decay configuration all 4-momenta are null vectors proportional to 
each other. 

3.2 Higher orders 

To formally compute the "exact" decay amplitude one has to replace in each tree diagram every 
boson propagator with the corresponding exact boson propagator, and each m-boson vertex 
with the corresponding m-boson proper vertex (i.e. one-particle- irreducible Green function), as 
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depicted in fig. (|^); there we have symbohzed each proper vertex by a blob. Diagrams involving 
ghost lines can be excluded for the same reasons as before. 




Figure 6: Exact amplitudes 

Using Property 2 it is easy to verify that if the external momenta are slightly off-shell, the 
momenta carried by the propagators in figg. (^), (|^) also are, and the scalar products of all 
momenta are of order uo] uo is the infrared regulator (with dimension of a mass) introduced in 
section 1. The exact propagators for massless particles in the infrared limit have to behave as 
the naive ones, i.e. they are of order cj"^. 

Let i?^, Ey., Eg and J-^, ly, Ig denote respectively the number of external and internal pho- 
ton,YM boson, graviton lines coming out of one of the diagrams in fig. (P). Let m^,m^,m^ 
denote the numbers of photons,YM bosons, gravitons coming out from the v^^ proper vertex 

appearing in the same diagram. Clearly, 



Ej 


= EK 

V 


-21^ 




Ey 


= EK 

V 


-2Iy 






V 


-2Ig. 


(16) 



Moreover, 

Np-Ip> e{Ep) (17) 

r if X = 

where 6{x) := < and Np denotes the number of proper vertices in the diagram where 

I 1 if X > 

at least one particle p {p being a YM boson and/or a graviton) comes out; this unequality 
follow from the fact that A'p = if and only if Ep = = Ip. 

The results of the next section (Property 10) can be summarized as follows, that 

= o(a;™7+4eK)-m;;+20(m;;)5o « ^ ^^g^^ 
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where in our notation o{uj^) will denote an infinitesimal or an infinite of at least order p in u, 
namely lim[o(c(j^)co'^^] is zero or finite. The overall tu-dependence of the diagram contribution 
D{u}) will be the product of the dependences of each vertex and each propagator: 



D{uj) 



Using equations (|T6|), the latter becomes 

D{uj) = o(^^--^-+4(^--^-)+2(^«-^«)), (20) 

where A^^ denotes the number of proper vertices in the diagram where at least one graviton and 
no YM boson come out. To estimate ^{Ny — ly) + 2(A^^ — Ig) let us distinguish two cases. If 
Ey = 0, then by colour conservation = for all vertices in the diagram, implying A^^ = A^^^; 
using formulae ioT p = y and p = g,we find A{Ny - ly) + 2{N'g - Ig) > 4e{Ey) + 2e{Eg). If 
Ey > 0, noting that (A^j, + A^^) = A'p, ly + Ig = Ip, where now p denotes either y or g, and using 
formulae (0) , we find 4:{Ny - ly) + 2{N'g - Ig) > 29{Ey) + 29{Ep) = 4^(Sj,). Summing up, 
4:{Ny — ly) + 2{N'g — Ig) > AO {E y) + 29 {E g) 5 Q ^ Thls expression depends only on the numbers of 
external bosons of the process, not on the particular diagram we are considering, therefore we 
find the following 

Property 7. - The amplitude T of a decay process with E^ external photons, Ey external YM 
boson and Eg gravitons satisfies the condition: 

rp ^ ^(^^E^-Ey+ie(Ey)+2e{E,)5f^ y ^^21) 

This formula is valid at any loop order in all particles different from the gravitons and at least 
at one loop order in the gravitons, because the matter action with a background metric is 



multiplicatively renormalizable [|1^, whereas at first order in the graviton loops pure QG is 
finite on- shell. 

Note that the RHS of formula (|21|): 1) is independent of the number of external gravitons, 
provided Ey > 0; 2) vanishes if Ey = 0. 

4 Ward identities 

In QED the proper n-photon vertices r^i "^"(p^, ■■■,p"') satisfy the Ward identity 

P..C^-^"(p^p^ ...,p")e,,(/)...e,„(p«) = 0, (22) 
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where is the momentum of the i-th photon and the corresponding polarization vector; 

this transversality condition amounts to the gauge invariance of any physical process involving 
n (incoming or outgoing) photons. 

In this section we first derive the identity above and its analogues for general Yang-Mills 
(YM) and Einstein (with A = 0) Quantum Gravity (QG) theories in the momentum config- 
uration of decay processes (compare with Property 2). Then we use them and a continuity 
argument to show that the proper vertex for any decay process with fixed external momenta 
vanishes in QED and QG, whereas it is finite in YM. The Ward identities are derived formally 
by using naive functional integration considerations based only on the gauge invariance of the 
classical action (not on its explicit form). In the case of QED,YM, their validity extends to 
the true (i.e. renormalized) theories at any order in the loops because renormalization pre- 
serves Ward identities. In the case of QG, their validity is guaranteed at any loop order in the 
matter fields and at least at one loop order in the gravitons, because the matter action with a 



background metric is multiplicatively renormalizable |13|, whereas at first order in the graviton 
loops pure QG is finite on-shell. 

We start by fixing the notation. Let Sinv{4>) denote the (local) action depending on the 
classical fields {cpi} and Ra{4>) corresponding (local) gauge generators: 

^^^^ = = 0, (23) 

<j<Pi 

We employ a condensed notation in which a capital index J is a collective index; it represents 
both discrete indices and a continuous space-time variables x. A repeated index implies sum- 
mation over discrete indices and integration over x. Explicitly, in the case of QED, YM, QG 
the fields 0/ include 

f A^{x), ip{x),%l){x) and/or ^{x), ^{x) in QED; 
y4^(x), + possibly i^^x), tp^lx) and/or v^*(x), ip'^{x) in YM; (24) 
hfj,u{x) + possibly any (pj considered in the two previous cases in QG; 

X G M'^ denotes the point in Minkowski spacetime, y4^(x), yl^(a;) the gauge potentials corre- 
sponding respectively to a f/(l) and a semisimple group G, ilj{x),ijj{x) (resp. (p{x), <^(x)) spinors 
(complex scalars), '?/'*(x), ^/'*(a;) (resp. ip'^{x),ip'^{x)) spinors (complex scalars) making up a finite 
multiplet belonging to some finite representation Rep{Cie{G)) (in the latter case (T°)* will de- 
note the matrix representation of the hermitean Lie algebra generators corresponding to A^), 
hf_iu{x) is the graviton field, ?7^jy denotes the Minkowski metric tensor (which plays the role of 
background metric) in cartesian coordinates, and g^u{x) = rj^y + Kh^^ is the the metric tensor. 
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The invariant actions Sinv read 

r -i V rf^x (F/^-^F^.) + S^at in QED; 

5... = I -i V rf'x (F'^ + Smat in YM; (25) 

[ Im^ d^x (A - Y^R) + Smat in QG, 

where F^^,, F^^ is the field strenght in QED,YM respectively, R is the Ricci scalar of the metric 
g^J.u, 9 '■= —det[g^^], /"'"^ are the structure constants of Cie{G) and e the coupling constant. 
Smat is the action of the matter minimally coupled to the gauge potential 

A^, y4°, h^y are respectively the gauge potentials for QED, YM, QG, with gauge transfor- 
mations 

5i,A^ = d^^ in QED; (26) 
5^A; = {D.O" ■■= d,C + er'^AlC in YM; (27) 

6^A; = A;d,e + edpA; inQG. (28) 

We omit for the sake of brevity the well-known gauge transformations of the other fields. 

The quantization of the theory (in a perturbative setting) is performed in the BRST for- 



malism |T2|, |TT[: the set of fields {(f)j} is enlarged to a set {$a} by the introduction of ghosts, 
antighosts and Stueckelberg fields, and we associate to the action Smv a gauge-fixed action S<i, 
depending on the gauge-fixing functional \E'. Index A, like /, represents both discrete indices 
and the continuous space-time variables x. Let Sgf '■= 5'^,($) — iSj„^(0); in QED and YM, Sgf 
can be constructed as Sgf = ■s\E', where s denotes the BRST transformation associated to the 



gauge transformations (|26| ) - (pSf). 

The generating functional Z{J) (depending on the external sources J) for the Green func- 
tions of the theory is defined by 

Z{J) := J P$ei[^*(*)+^-'*^], (29) 

where P$ is a gauge invariant functional measure, transforms under diffeomorphisms as 
the appropriate tensor density. 

Strictly speaking, in the case of QG an action Smat containing a spinor contribution requires the introduction 
of vierbeins as dynamical variables instead of the metric. However, the considerations of this section hold also 
in that case, since they are based on the gauge tranformations (^8|) of the metric, which can be obtained from 
the gauge transformations of the vierbeins. 
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By performing a gauge transformation — + S^(f) of the dummy integration variables 
in the RHS of eq. ( ]29|) the integral Z{J) remains the same (the Jacobian is 1), implying the 
Ward identities 

= 5^Z{J) ='-J P<I>[J^55$A + 55S'GF]e^[^*(*)+^''*^l, (30) 
or, in terms of the generating functional W{J) := ^ln[Z{J)] of the connected Green functions. 







W{J) + disconnected terms. 



(31) 



The disconnected terms are absent when evaluating the Green function on any decay process, 
since in this case only one initial particle is present. Therefore, as far as we are concerned. 







W{J). 



(32) 



In order to obtain the Ward identities for the proper vertex functions we introduce the usual 
Legendre transform r(<i>) := [W^(J) — J^^a]\j=j(^), where the function J = J($) is obtained 
by inverting the relations $^ = the new independent variables are the "classical fields" 
Consequently J^(l>) = 



From identity 
proper vertices F 



we draw the following Ward identities for the generating functional of 
(5F 







5$, 



6^^A + 555gf($) 



(33) 



Actually, we are interested in the Ward identities for the proper vertices having only physical 
gauge bosons as external (incoming or outcoming) particles. The physicality condition is best 
imposed in momentum space. The proper vertex rl^'''"'{x^,x'^, ...,x^) with n external gauge 
bosons 6i(a;*) (in configuration space) is obtained from F through differentiation, 

5"F 



^12...n/l 2 



(a;\a;^...,a;") 



6bi{x^)...6bnix' 
where we have introduced the short-hand notation 

4, inQED 
i^, in YM 



(34) 



'I>=0 




1,2, ...,n. 



(35) 



Vi', or in QG; 

The RHS has automatically the required boson symmetry in the identical particles, e.g. if all 
the fej's are the same type of fields 



,x 



F 



12...n/l 



(x , X , X ), 



(36) 



■'-■'-Alternatively, one could perform a BRST transformation; the resulting Ward identities would be the same. 
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where {ii,i2, ...in) is a permutation of (1, 2, n). On account of the translation invariance 
,...,x"') = r^" "(x^ + a, + a), its multiple Fourier transform can be written as 
ri'""(p^, ...,p"')S^{J2 P*); it contains a Dirac-5 implementing the total momentum conservation. 



i=l 



Here and below our conventions for the Fourier transform will be f{p) := / ^ ^p-xj!^ 
f{x) = J d'^p e^^'^ f (p) . As a consequence of the general relation 



d^x 



-ip-x 



.. 6F 



one finds 



(2vr) 



-An 



6F 



5"r 



\i=l 



6hi{~p^)...6hn{-p" 



(37) 



(38) 



<I>=0 



Differentiating relation (|33D with respect to 6i(— p^), ...,bn{—p"') and setting thereafter $ 



0, we obtain 



d^q 



'o(g) 



/l=l 



+ 



56i(-pi)...56„(-p") 



(39) 



$=0 



In fact, only the terms with $^ = 6 in the first term in eq. (|33| ) contribute to eq (|39|), since 
when <I>A 7^ 6 then ^^^^^^^^ | |,=o = (indeed, for any $a S^^a is of degree > 1 in $^). 

To get identities involving proper vertices with physical external bosons we will have to 
contract their Lorentz indices with the ones of transverse polarization tensors/vectors (we will 
choose them with well-defined helicity) e^(p^)...e"(p"), where 



e{p) = e^{p) :-- 




when 6 = A^, 



when b = h. 



with e^{p)p'^ = 0. 
Now it is easy to realize that in all cases the following property holds: 



et{p)et{p)) 



(40) 



6bi{p^)...6bn{p'' 



i{-p')...e:{-pn = 0; 



(41) 



$=0 



where contraction of the Lorentz indices hidden in the symbols l,...,n and e^, ...,e" is under- 
stood. In fact, the terms of non-zero degree in the ghosts contained in 6^Sgf vanish after 
setting $' = 0; the other terms depend on the longitudinal modes of the bosons, and vanish 
after contraction with the polarization vectors/tensors. We prove explicitly this statement in 
the appendix, for the Feynman (harmonic) gauge fixings. 
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Introducing the notation 

r^-"'-" := r^-'-" ■ e\ (42) 

where again contraction of the Lorentz indices hidden in the symbols i and e* is understood, 
the Ward identities (143|) will therefore reduce to 



(43) 
l>=o 



h=i \ j=i, 

The identity above is one essential ingredient that we need in order to prove the main 
property of this section. In order to formulate this property, we need now a notion of "vicinity" 
of a "decay configuration" parametrized by one regularization parameter uj. Therefore, we 
introduce some useful definitions. 

A configuration cj-converging to the decay configuration (^*, £*)i=o,...,n {'^ > 0) is a one- 
parameter family (A;*(u;), £:*(co'))j=o,...,n such that {uj) ■ k"^ {uo) = 0, k\uj) — k''- = o{uj), e''-{uj)—i'^ = 
o{iu), k^ ■ k^' = o{uj'^) Vz, z' = 0, 1, ...,n. Examples of these families will be given in formulae 

(P, m- 

It is easy to show that in the mentioned hypotheses the 3-momenta are in general no more 
coUinear, but form angles <, u; consequently, 

{either o(l) 
e\k') ■ y = o{uj). (44) 
or o{uj) 

We are now able to prove the following fundamental property of the vertices, which is the 
main result of this Section and adds to the kinematical properties of Section 1: 

Property 10. - On any configuration (A;*(c<j), e*(u;))j=o,...,n o^-converging to the decay configu- 
ration (A;\£*)i=o,...,n 

Ti^^'\k\ r) = o(a;"+^) in QED; (45) 
peOao...e"a„^^o^ fc") = o(cu^-"'^) in YM; (46) 
rf+f " (A;°, A;") = o(u;™^+^(™i')(4-™y)+2e(m.)C ) in QG. (47) 

where in the third equation m^,my,mg denote the number of external photons, YM bosons 



and gravitons respectively (m^ + niy + nig = n + 1), and 9{x) 



if X = 

1 if a; > O' 
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Proof. 

The claim is evidently true when n = 0. In fact, Fj*" oc {k^Y° in QED, YM, but this 
vanishes since momentum conservation imposes the condition = 0; in QG still it could be 
pAto^^o _ ^Qjigf X ryW^'o^ 1^^^ |^]-^jg vanishes after contraction with e^"'^" (which is a traceless tensor). 

The rest of the proof is by induction and divided in three parts. Let us assume that the 
claim is true when n = m — 1. We will prove that it is true when n = m. For the sake of 
simplicity, we explicitly prove the claim (|47|) , which is the most general possible, in the simpler 
case = = niy, 

at the end of this section we will briefly sketch how the proof goes in the general case. 
Part 1 Here we prove the equations 

peO...»-,..e»+^..e"^^o^ k^k^^ = in QED; (49) 

^e^ao...e^-^a,.,,^,a,,e^+^a.+,...e-a^ YM; (50) 

''^'"'^"'■■^"(^°, k^k;^ = o{u') in QG. (51) 

We drop in the sequel the tilde and write A^j^, A^^, g^^ instead of ^4^,71°, g^^. We treat separately 
the cases of QED, YM and QG. 

- QED. From d^^A^ij)) = ip^Hp) (eq. (p6D ), and eq. (^3|) , from differentiating w.r.t. q it imme- 
diately follows 

pJoC+t-'^7/,p\-,P")=0 (52) 

n 

(we have factored out 6^{Y1, P*))? whence formula (|49|) follows at once (using boson symmetry), 

2 = 

if we choose so that the sets {p^, p"}, {A;°, A;"} coincide. Actually we can derive directly 
from eq. (|39|) the stronger property 

A:; r^--^-'^"(fcO, k\ r) = 0, n>2 (53) 

- YM. From 

= ^P,C{P) + er"' J d'q A^p - q)C{q) (54) 

(eq. (|2^) in momentum space), and from differentiating formula (|43| ) (with n = m) w.r.t. $,{p^), 
it immediately follows 

'irP ■pMoao,eiai,...,e'"a„/ 1 m\ , 

'P/iQ-*- m+1 \P 1 P 1 P J "I" 

m 

^Y^^fhaia^Y<'^au...e'-^ai_^,e%,e^+^a,+^,...,e^a^ = (55) 
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m 



(again, we have factored out P*))- This formula holds for any configuration I] = 0, 

4=0 i=0 

e*(p*) -p* = 0. On a configuration cu-converging to the decay configuration we deduce from the 
induction hypothesis that the second term is o{uj'^~"^). 

- QG. The gauge transformation (pSj) in momentum space reads 

Si;9fw{p) = i J d^r{gp^{p - r)r^^P{r) + gp^{p - r)r^^P{r) + ^^(p - r)rp5(^.^(r)}, (56) 
implying 



Moreover, we note that 



Sgapjp) 



5\p + q) « + « . (5^ 



After differentiation w.r.t. ^'^"(p^), Eq. (|43|) with n = m reads: 

= cT^•••''^'"(p^p^•••,p")2(/).^..o 

m 

+ E [^:n■••^"'''^'^'^"'■■'="(•••,p'-^p° ••04(p\(e'^).(^'^) 



h=l " 

m , m , 

(once again, we have factored out 5^( P*))- This formula holds for any configuration = 0, 

1=0 i=0 

e*(p*) -p* = 0. On a configuration tu-converging to the decay configuration we deduce from the 
induction hypotheses ([51]), ( P8| ) that the second, third terms are o{uj'^), which proves eq. ( pT] ) 
for n = m. 



Part 2: We prove the factorization formulae 

rf+i'"(^°,-,^") = E^Mi...in^'"*^-^'"-^'" in QED, (n+1) even: (60) 
p 

peOao...e"a„^^o^ F) = ^ A^^^X^'"'' + in YM, (u+l) even;(61) 

p 

r:V"(^°, A:") = E A.on.....„..„+i^^°"-^"""""^ + o(^) in QG. (62) 
p 

and 

n 

r:V-"(fc°, r) = 5: ^ ■ e^")i?^^^^..i?'"-^^" in QED, (n+1) odd; (63) 

P jo=0 

n 

peOao...e"a„^^o^ P) = ^ ^ ^oT.'.in " l'^^" " s'"') E'''\ . .E'-''" + o{uj^-'') in YM, (u+l) odd; 

P jo=0 

(64) 
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where: 

1) J2 means the sum over all the permutations P [P = is a permutation of 
p 

(0, 1, n) in QED and YM, whereas P = {io, ii, Z2n+i) is a permutation of (0, 1, 2n + 1) 
in QG); 

2) the A's are scalar functions depending on the scalar products ■ (and, in the Y.M. 
case, on 2m Lie algebra indices Oj); 

3) we have introduced the shorthand notation 

E'^ := (e^ ■ k' ■ F - e' ■ F ■ F) . (65) 

In the RHS of eq. {^) it is tacitly understood that e^^+i = e"^', k^'+^ = k"^', s = 0, 2n. 

We prove explicitly the first three (the proof of formulae (|63D, (|6^, is completely analogous): 
let + 1 = 2m. We look for the most general r(^J^^"^"^^(A;^, k"'^^) satisfying: 

1) the constraint 

rir'''^-^"'-^"(A:°,...,r)fc;^ = o(a;'^) (66) 

in any configuration (/c*(co'), £:*(u;))j=o,...,n ci^- converging to the decay configuration (A;*, e*)i=o,...,n; 

2) symmetry under any replacement (/ij, k^) ^ (///, fc'), i,l = 0, n. 

If we set o{uj'^) = this amounts to solving eq. (^9|) equipped with boson symmetry for 
the (n+l)-photons vertex function of Q.E.D.; if we set c? = 3 — n, this amounts to solving eq. 
( |50D equipped with boson symmetry for the (n + l)-gluons vertex function of Y.M., provided 
we understand an implicit dependence of r„+i on the Lie algebra indices Oj and remind that 
the latter have to be permuted along with the indices /ij and the momenta /c* when boson 
symmetry is imposed; if we choose n + l = 4r, d = 2, and add the additional symmetry 
conditions = k'^\ 6^*+^ = e^* (i = 0, ...,2r — 1), this will amount to solving eq. (^Tj) 

equipped with boson symmetry for the 2r-gravitons vertex function of Q.G. In this way, we 
can formally deal with eq.'s (|49|) , (|50|) , (|5T|) simultaneously, by just dealing with one. 

The dependence of r^5^']^''^"(A;°, fc") on Lorentz indices can only occur through the metric 
tensors r/'^'^^j and the 4- vectors k^^'. Compactly, the most general dependence can be written in 
the following way 

pMo..^..^„ _ ^ ^0 ^^^^^ +Y,B^kk r]...rj +... + J2B"' , (67) 

m times (m-1) times ^"^ 

where the S's denote Lorentz scalar functions. For our purposes, it will be more convenient 
to expand F^+i in terms of the 4-vectors k^' and of the tensors E^'^^{k\ y) := ri'^^^^^k^ ■ k^ — 
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[k'^Y^ {yY^)^ which satisfy the relation 

= = (F)^^E^'^^ (68) 
The general expansion (|^) can be replaced by 

m n 
P 1=0 jo, -121-1=0 

(69) 

where J2 nieans the sum over all the permutations P = (zq, h, ■■■,in) of (0, 1, n) and A^-J]^^-'^' 
are scalar functions depending on the scalar products k"^ ■ k^ (and, in the Y.M. case, on 2m Lie 
algebra indices a^). 

We have introduced a quite redundant set of scalars {A^- J^f'-'^'} to make formula (|69D more 



compact. The set is redundant in the sense that A i^",-^''"' and A?' ^°r^''-' will both contribute 
to the same term [yoyra ___Q^j2i-iY^2i-i E^''ii^''ii+'^ ...E^^"-"^^^" in the expansion (||), whenever 

1) there exists a permutation P21 of 21 objects such that (io, "^i, ■■■,i2i-i) = AK^O; "^i; '^2«-i), 
(jo, jl, j2Z-l) = P2i{jo,ji, j2«-i); 

2) («2h^2/+i, •••,»in) = -Pn+i-2/(^2/,^2i+i, •••,^n), whcre P2m-2« IS a pcrmutatiou of n + 1 - 2/ = 
2m — 21 objects which is the product: 2. a) of transpositions between the (2s)*^ and the (2s + l)*^ 
object (s = 1, m — /); 2.b) of transpositions between different pairs (2s, 2s + 1), (2r, 2r + 1), 
r, s = 1, m — /. 

We are free to set A^^T^''-' = A'' ^°r^''-' in these cases. 

*0---Jn to ■■■In 

Finally, boson symmetry ( p6|) implies that the scalars A' satisfy the relations 

{j ii h = i 
% ii h = j (70) 
h a h ^ i,j 

for any pair of indices 

Plugging the general expansion (H) into Eq. (H) and using relation (|68|) we find 



m n 



X 



P' 1=1 jo,-j2l-l=0 

■ fc^o) e'^ ■ k^\..€'^'-^ ■ y^'~^E'^''^'+\..E'^-^'", (71) 

where means the sum over all the permutations P' = (zi, in) of (0, 1, z — 1, i + 1, n), 
whereas 



m n 



1'; io--i2i-\ I 4^; 30—321-1 I I /i'; io---j2i-i 

^jl...i2i-lM2i---«n 



pe"...e"/7,0 _ VV V A^', 30-32l-i , 30-32l-l , , 

P' 1=0 30,-321-1=0 

(e* ■ F'o) ■ . k^2i-iE'^2ii2i+i^^^E''^-iin_ (72) 
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Note that the term / = has completely disappeared from the sum in eq. ([7I|), due to eq. (|6^). 

Let us fix the xyz axes so that = (/cq, 0, 0, k^) [according to property 2 this implies 
k^ = A-^(fco, 0, 0, fco), j = 1,2, ...ra]; we can always assume that the polarization vectors are 
real and have the form = (0, cos6^, sen6^, 0). We now start exploiting the available freedom in 
the choice (1) of the angles 6^ characterizing the polarization vectors e*; (2) of the configuration 
C(j-converging to {k\i^)i=o,...,n- A family of possible choices of the latter is 



k' = k' + uUe' 



£> ■= (0, -sene\ cos9\ 0) 



i = 0, 1, ...,n. 



(73) 



the family is parametrized by the 2n + 2 parameters (6*, 6'*), which are only constrained by the 

n . „ j n , n ^ , 

condition Yl, b^^' = (so that J2 = = 0). As a consequence 

i=0 i=0 i=0 

. y = ~uj%'Vcos{e' - e^), e' ■ y = -ub>sin{e' - e^) e' ■ = -cos{e' - e^) (74) 

(e' ■ k' ■ F - ■ F ■ k') = uW (75) 
By plugging these {k\e^) into Eq. (^) we find 



m n 



30---j2l~l _|_ j^l; jo---j2l-l I \ A^' i0---i2i-l 



1 I -L A ' 

111. ..In ^ ■■ ^ ^Zl...l2l-lll2l...in 



EE E E[4; 

P' 1=1 i0,-j2i-l=0 i=l 



(76) 



The coefficients in the square brackets can depend on the angles 6^ only through the cosines 
cos{6^ — 6^) (since /c* ■ y = —uj'^UV cos{9'^ — 0-'))] since the above equation has to hold for all 
6'*'s then all terms in the square brackets have to satisfy the relation 



30. ..321-1 _|_ /I 



l\ 30---321- 



1 _l_ _|_ 4'' io— 321-1 

-t- ... -t- ^i^,„i2i_iii2l...in 



d-n-2\ 



Z = 1, m 



(77) 



independently. 

Replacing the above results in formula (^) we find the factorization formula 



p 



whence formulae (|60D, (|6TD , ( |62D follow. 



(78) 



Part 3: On any configuration (A;*(co'), e*(c(j))j=o,...,n ix'-converging to the decay configuration 



Ck\i') 



i=0,...,n 



we have E^^ = o(c(j^). To prove formulae 



61) , (HHD it remains to show that the 
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scalar functions A's appearing in eq.'s (|60|), (pT|), ( |62D can show poles in oj at most of degree so 
high to yield the global cj-dependence reported in the former formulae. For this purpose we use 
a continuity argument, i.e. we argue that the claimed tu-dependence is the only one compatible 
with equations (§3), (|6TD, ( |62D if we require the LHS to be independent of the particular 
configuration £:*(u;))j=o,...,n w-converging to the decay configuration (A;*, e*)j=o,...,n- 

For the sake of brevity we continue to use the factorization formula (|78D to deal at once with 
all three cases. We choose two different multi-parameter families (A;*(u;), {k^{ijj),e^{uj)) 
of configurations cj-converging to the decay configuration, and we require that 

eO(.o).....(..) ^ linirfif)--(^")(P,...,r) (79) 



In the xyz axes as before, the first is the family ([73D, the second is 

F + u;(0,0,c\0) t ■=[{kiif + {ue)X^{<d,Q,kl,ue) i = 0,l,...,n, (80) 

where = 0. This implies in particular fc* ■ y = —cu'^dc^. 

i=0 

With the first family we find 

r^^-"(A;°, fc") = c."-*-i6°...6" Y: ^n...^. + oiu;'-'). (81) 

p 

Now we specialize our discussion to the case of QED and QG, where d — 1 > 1, so that 
the second term vanishes when — > 0. Let us consider per absurdum the hypothesis that the 
functions A's have poles of degree (n+l) in u. In order that the RHS has a limit independent 
of the 6"s when uj 0, the A's must have the form 



^io«i...«„ 



L P 



-1 

(82) 



where aiQi^...j„ are constants, so that 

r n -1 

X const. (83) 



A"' 



cj"+l6n...6" 



On the other hand, plugging the family ( pOD into eq. (^2|) and replacing the result into 
formula (|78]), we find 

rf <r. ■-•".(P. .... n - c„... .^(li,^"-AJ^-,^-A, (84) 



p 



where we have defined d^ := This expression depends on the choice of the coefficients c*, i.e. 
depends on the way the family (k^u) , e^u)) approaches (k'' (u) , e\uj)) , against the hypothesis. 
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In a similar way, one can exclude the hypothesis that the functions A's have poles in uj of degree 
> (ri + 1), otherwise the RHS would diverge to either +00 or —00 according to the way the 
families approach the decay configuration 

Summing up, we have discarded the possibility that the A's have poles in lo of degree > n+1, 
so that consequently in QED,QG 

r^^-f (A;°,...,r)=oM (85) 

In QED we can improve the bound ( ^5|) into the stronger bound (|45|) . In fact, if one plugs 
the general expansion (^) into eq. ( p^D [instead of eq. (pB|)] and argues as in part 2, one ends 
up with a stronger form of the factorization, 

C+r^"(fc°,...,r) = Y.M,i,,„,^E^^^^^^...E^^--^^^-. (86) 

p 

Looking at the Feynman diagrams contributing to each order in the loops to Y^^^"^ {k^ ^ fc"), it 
is easily realized that they are continuous and finite for all values of /c*'s, since the fermion/scalar 
masses are infrared cutoffs [see fig. (^]. Hence, the scalars A cannot have poles in k^-k\ because 
otherwise at least the terms Aj(,ji.,.j^(/c*o)''*i ...(F")'^'"-i (/c*"~^)^''" would diverge. The A's have 
dimension [mass]^"^*-"^^-', since r„+i has dimension [mass]^"*-"^^''. This can be accounted for 
without introducing poles in ■ y, but using the mass parameters of the charged particle 
interacting with the photon. For instance, if the only charged particle is a fermion with mass 
m, then A = m^~^^"+^)o(l). We have completed the proof of the claim (|45|) . 

In QG the o{uj) in the RHS of (|85D can be improved into a o{uj^), since T'^^f"{k^, fc") can 
be only of even degree in cj, if we assume that the proper vertices depend analitically on the 
momenta /c*. This follows from formula (|4^), because the LHS of eq. ( ^5)) has to be a function 
of the Lorentz scalars /c* ■ k\ e^{k^) ■ y , of even degree in the latter. This completes the proof 
of the claim (^). 

In YM formula (|7S|) and the continuity argument do not exclude that there exists a limit 
lim r^V{^"(/i;°, fc") =: L 7^ independent of the way the family {k^ (u) , e^u)) approaches 
[k"^ {uj) , {uj)) . In fact, if the functions A's have a pole of degree > {n+1) in uj, the second term 
in formula (^) (which in principle can be finite or divergent) could compete with the first, 
and r^°_^^[^"(fc°, fc") could have a family-independent limit even though the first term has not. 
This is exactly what happens with the 4-gluon proper vertex, as one can already check at the 
tree level 
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By an explicit analysis of the general expansion (p9D one can easily realize that a family- 
independent limit L G R U {±00} can be obtained only if equation ( ^61) is satisfied. 

Finally, the proof of the general claim (^) can be done by an induction procedure in the 
number of external photons (resp. of YM bosons) which mimics the one sketched so far for 
QED (resp. YM), with the only difference that as starting input we do not use the value of 
proper vertex with zero photons, zero YM bosons and zero gravitons, but the proper vertex 
with nig > gravitons or rriy > YM bosons (resp. with rUg > gravitons or > photons). 

We have thus completed the proof of property 10 <C>. 

5 Concluding remarks. 

We have seen that the decay probabilities for the photon, the graviton and the Yang-Mills boson 
all vanish (perturbatively). The decay amplitudes involving only photons and/or gravitons are 
themselves zero; we have first shown these properties by a simple power counting argument 
and then proved them rigorously through the Ward identities, assuming only continuity of the 
Greens functions in the infrared limit. In the case of the Yang-Mills boson, the amplitude does 
not vanish in the infrared limit (more precisely, it diverges if m > 5 out of n+1 external particles 
are YM bosons); the decay probability is however suppressed by the phase-space factor. The 
latter is the only case in which we have needed an infrared regulator. 

In this final Section we would like to comment on the relation between our work and the 
classical literature on infrared divergences in quantum field theory. 

For the reasons just mentioned, even in YM theories we do not need to average (d la 
Block and Nordsieck [0]) over sets of states degenerate in the energy, like in the Kinoshita-Lee- 
Nauenberg theorem , in order to build finite physical transition probabilities out of divergent 
amplitudes. 

However, this might be necessary for other theories, not explicitly considered here, where 
the divergences of the amplitudes are sufficiently bad. In the latter case it would be nevertheless 
important to keep in mind some peculiarities of the decay of massless particles compared to 
what one usually finds in the literature 0. The physical processes explicitly considered in the 
literature are either scatterings, or decays in which the initial particle is massive. The case of 
a decay process where all particles (including the initial one) are massless is not considered. 
In a theory including massless particles, the Green functions may diverge if (1) some of the 
external particles are soft/collinear massless ones and/or (2) if massless particles appear among 
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the internal ones occurring in the corresponding Feynman diagrams (e.g. in loops). The study 
of these divergences is usually performed by attributing a small mass m to each kind of massless 
particle in the theory and then studying the limit in which m goes to zero (as already recalled, 
their elimination from "physical" transition probabilities is obtained by building up the initial 
and/or final states as a mixture of degenerate states of the energy before performing the limit; 
see Ref.'s g). 

On the other hand we note that, while in the scattering processes or in the decay processes of 
a massive particle the coUinearity of some massless external particles is one of the kinematically 
allowed configuration [so that it makes sense to study the divergences of the Green functions 
in the limit when these external momenta become coUinear while remaining on-shell (null)], in 
the decay of a massless particle the only allowed kinematical configuration is that in which all 
the external particles are collinear. Therefore, the divergences, if they occur, characterize all 
the kinematically allowed configurations. Moreover, the latter tipically appear already at the 
tree-level (consider e.g. T5 for the YM theory), as one can immediately check by summing all 
relevant tree diagrams. As for the IR regulator, the one based on the attribution of a small 
mass to the external particles is unsuitable because it forbids the decay of a particle into other 
ones of the same kind (as we have already noted in section L3). Our regulator, based on a 
small-frequency external source, bypasses this difficulty while having the nice feature of being 
physically intuitive. 

As mentioned at the end of Section 2, a partial decay probability r„ different from zero 
can be only obtained when the square amplitude is proportional to a sufficiently high negative 
power ofuj. If we admit (as is generally true in perturbation theory) that the coupling constants 
appear in the numerator, this means that the amplitude must contain a coupling constant with 
positive mass dimension. 

One of the few theories we are aware of, in which such a coupling occurs (besides the X(j)^ 
theory; compare Section 1) is gravity in the presence of a cosmological constant. In this case 
the action of the gravitational field is written as 

S = J d'x ^) [A - R{x)] (88) 

or, redefining the metric in the form g^u{x) = rj^y + Kh^^{x), with k = vTGttG, 

S= I d^x ^1 + Kh + k'^K^ + K^lfi + ... [4 - R^^Hx) + ...1 (89) 
J J 

We have denoted symbolically with h, h"^, ... in the square root terms which are linear, 
quadratic, cubic ... in h, omitting the indices and the exact algebraic structure. RP'\x) denotes 
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the part of the curvature quadratic in h. The term n^h^, when is multiphed by A/k^, gives rise 
to a vertex nAh^ which couples three gravitons with a couphng constant kA of mass dimension 
1 (unhke the corresponding three-vertex of the pure Einstein action, which is proportional to 

and contains 4 four- momenta, so that the infrared processes are strongly suppressed). 

Although the decay amplitudes involving this new three-vertex are suppressed at the tree- 
level because of helicity conservation (Property 3), it can be used to construct gravitonic loops 
with n external legs. The amplitudes will be proportional to positive powers of kA and - in 
our regularization scheme - to negative powers of u. This means that r„ would be finite in 
the limit uj 0, or even diverge. But we should not forget the terms which are linear and 
quadratic in h in the square root of eq. (|89|). In particular, the quadratic term gives rise to 
a graviton mass (if A < 0) or to instability (if A > 0) 0|. In the first case, we end up with 
gravitons which are not massless any more, so that all our preceding formalism does not apply. 

It is known that the cosmological constant A, although possibly very large in principle, is 
limited by astronomical observations to be less than |A| < lO^^'^G"^ (in order to explain this 
vanishing, many mechanisms have been proposed Therefore it seems that the idea of a 

decay induced by the presence of a cosmological constant can be excluded on the basis of the 
empirical evidence. 



However, in the non-perturbative quantum Regge calculus |T5| the effective value of the 
adimensional product \A\G depends on the length scale and vanishes according to a power law 
as the energy scale fi goes to zero: 

|A|G ~ (M^ (90) 

If we admit that the average lattice spacing Iq is of the order of the Planck length , then the 
constant A can be non-vanishing on small scales, leaving the graviton massless at large scales. 
This might change the situation, but clearly at the present stage of knowledge these are still 
speculative hypotheses. 

We wish to thank M. Abud, P. Van Nieuwenhuizen and G. A. Vilkovisky for useful discus- 
sions, as well as the A. Von Humboldt foundation for financial support. We are grateful to D. 
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Appendix 



We prove eq. (|4T|). In the case of QED with (for instance) Feynman's gauge-fixing J d^x^d^A^)"^, 
the LHS is zero when n > 1 because the gauge variation of the gauge-fixing above is of first 
degree in A^, and is zero in the case n = 1 because 

- -PiiPimpi), (91) 

a 

vanishes after contraction with the polarization vector £^(pi). In the case of YM with (for 
instance) Feynman's gauge-fixing ^ / d^x{d'^Afj,)'^, the LHS is zero if n > 2 because the gauge 
variation of the gauge-fixing above is of second degree in A^; if n = 1 it is zero for the same 
reason as in the preceding case (|9ll) ; if n = 2 it is zero because 

2 



a 



-{pinP2rc{pi+P2)r'''', (92) 



vanishes after contractions with the polarization vector (pi)£^^2(p2)- In the case of QG with 
harmonic gauge-fixing J d^x{d^h^jjY we have 



5, 



^ j d'x{d^h,,f'\ =~J d'p{p'^h,,{p))f{(^ + &){-p) (93) 



(here " means Fourier transform). When some .(-p ) ^'^^^ P'^h^,y(p) we get a factor 5^{p + 
Pi)[^il'K'' + ^u'^'il'] (see formula (|58|)), which gives zero after contraction with the polarization 
tensor e^^'^^{pi). 
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